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The superconducting state around a well-isolated nonmagnetic impurity in the high- Tc su- 
perconductors is studied using the two-dimensional t-J model. The spatial dependence of the 
order parameter and the local density of states are obtained from the numerical diagonalization 
of the Bogoliubov-de Gennes equation derived using the Gutzwiller approximation. We find a 
zero-energy peak in the local density of states around the impurity. Different from the previous 
results on a vortex or surfaces in the t-J model, the splitting of the zero-energy peak is not 
found. The zero-energy states corresponding to the zero-energy peak is approximately localized 
around the impurity. 
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impurity, zero-energy peak 



The role of impurities in superconductors has been 
studied theoretically since the establishment of BCS 
theory. In the conventional (s-wave) superconductors, 
nonmagnetic impurities have only little effect on the 
transition temperature, as understood from Anderson's 
theorem. EP On the other hand, they have strong ef- 
fects on the superconducting properties of unconven- 
tional (anisotropic) superconductors in various ways de- 
pending on the anisotropy of the pairing state. There 
have been many studies on the nonmag|iietic impurity 
effects in connection-Yfithnheaygz^f^imonBu'cl' and high- 
Tc superconductors.! B'QB&Byyy 

Regarding high- Tc superconductors, a large number 
of experimental results and theories support the d^2_y2- 
wave pairing state. The significant difference between 
the dx2_y2--wsive pairing state and the conventional one 
is that, in the d2:2_j^2-wave state, the order parameter 
(pair potential) changes its sign with 7r/2 rotation. Thus, 
the nodes of the energy gap exist on the Fermi surface 
along the kx = ±fcj, directions in the Brillouin zone. The 
existence of the nodes implies that there are residual 
quasiparticles around the nodes. Thus, tha|CL2_j^2-wave 
pairing state has an intrinsic instability .c3tB This in- 
stability could show up in nonuniform systems, such as 
those with vortices, surfaces or impurities. In these cases, 
we expect an interference effect of the quasiparticles be- 
cause the quasiparticles feel the sign change of the pair 
potential when they are reflected. Actually this inter- 
ference effect causes the zero-energy states, which can 
be detected as the zero-energy peak (ZEP) in the-scan- 
ning tunneling spectroscopy near the vortex corala' and 
surfaces Jl3) 

In this paper we study this interference effect in the 
nonmagnetic impurity problem. To date, there have been 



some works on the nonmagnetic impurity in the d^^-y^- 
wave superconductivity. It has been shown that the low- 
est eigenvalue approaches zem^ the impurity scattering 
approaches the unitary limit .&E2P This is the zero-energy 
state which is located around the impurity. In this pa- 
per, we show for the first time that this state actually 
gives the ZEP, which can be observed in scanning tunnel- 
ing spectroscopy around impurities. In order to obtain 
the local density of states around a well-isolated impu- 
rity, we study fairly large systems. Furthermore, we will 
show that the zero-energy state is localized around the 
impurity. Our result gives an answer to the inconsistency 
between the previous theories, that is, whether the local- 
ized staie_has a slosv-decay {^ l/r) along the nodes of 
the gapB'lB) or notMEJ 

Finally, in contrast to the vortex coresH and sur- 
faces ,113' we show that the extended s-wave component 
is not induced around the impurity, and thus there is no 
splitting of the ZEP. 

As a model for high-Tc superconductors, we use the 
two-dimensional t-J model. Quasiparticle states near a 
nonmagnetic impurity in the dx2-y2-'wave superconduc- 
tivity is studied for the case of a unitary scattering limit, 
which is discussed in the previous works and is relevant 
to the Zn-doped case. The spatial variation of the su- 
perconducting order parameter around the impurity is 
determined self-consistently. 

The Hamiltonian of the t-J model with a nonmagnetic 
impurity is written as 
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where {i,j) means the summation over nearest- neighbor 
pairs, and the impurity potentials V/™'' are nonzero only 
on the impurity sites. Hereafter, we use t as an en- 
ergy unit. Since the constraint of no double occupancy 
is imposed on this Hamiltonian, it is difficult to carry 
out analytical calculations even in mean- field theories. 
Here we consider T — Q variational theory. The varia- 
tional wave function is PG|BCS(Aij)), where Pq is the 
Gutzwiller projection excluding the double occupancy, 
Pg = H,(l - n^^Uii), and |BCS(Aij)) is the BCS mean- 
field solution with the site-dependent order parameter 
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evalute the variational energy with the GuJ^wilk 



In order to 
ler projec- 
tion, we use a Gutzwiller approximation,lfd'EiP in which 
the constraint is taken into account as a statistical av- 
erage. This approximation has been extensively studied 
in connection with the variational Monte Carlo (VMC) 
calculationEB' For the two-dimensional t-J model, it 
has been shown that the Gutzwiller approximation and 
the iZMC calculation give very similar variational ener- 
giesJfil'' 

Using the Gutzwiller and a mean-field approximation, 
we obtain a Bogoliubov-de Gennes equation and a set 
of self-consistent emrations in a similar way to the BCS 
mean-field theory JlS-' 
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where r runs as vectors pointing to the nearest-neighbor 
sites and gt, Qs are the renormalization factors in the 
Gutzwiller approximation given as functions of the dop- 
ing rate 5 = 1 — n: 

25 4 



9t 



9s 



1 + 5' "" (1 + ^)2- 
Self-consistent equations are 

= -\ E(""*«; + ""<*) tanh ^^ 
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= -- Y^iurW} - vfvn tanh ^, 



(5) 



In the following, we assume S,ij] — S,iji = S.ij and that 
Ay is a singlet pairing, i.e., Ay — Aji. Since we consider 
the dilute limit of the impurity concentration, /^ is fixed 
to the bulk value /io determined without impurities. 

In the present calculation, we regard the Nl x Nl 
square lattice as a unit cell of which the impurity is lo- 



cated at the center {Nl is odd. See Fig. 1). We use 
periodic boundary conditions in x- and y-directions and 
assume a translational symmetry of Ay with respect to 
this unit cell. Then we can make use of the Fourier trans- 
form of the Bogoliubov-de Gennes equation. In order to 
study the effects of a single impurity, it is necessary to 
take the size of the unit cell large enough. Specifically, 
Nl should be chosen so that the typical V-shape local 
density of states is obtained on the sites at the edge of 
the unit cell. For the Fourier transform, we take the 
number of the unit cells Nc — 72 (8x9). 




Fig. 1. The unit cell for the impurity problem (enclosed by the 
double lines). Ni^ is odd. The black circle represents the im- 
purity site. A, B and C represent sites where the local DOS is 
calculated. 



We solve numerically the Bogoliubov-de Gennes equa- 
tion and carry out an iteration until the self-consistent 
(4) equations for Ay and ^y are satisfied. We take J/t = 
0.2, J^/t = and the doping rate 5 = 0.05 throughout in 
this letter. These parameters are approximately consis- 
tent with realistic ones obtained by Hybertsen et al.c3 
The impurity potential is taken V™^^ = lOOOt on the 
impurity site so that the impurity scattering is in the 
unitary limit. 

Figure 2 shows the obtained order parameters in the 
unit ceU for Nl = 27 and 5 = 0.05. We define Af and Af 
for each site by averaging A^ ^^^.'s on the two opposite 
bonds of X- and y-directions. Furthermore, the obtained 
Af and Aj' are decomposed into extended s-wave and 
d-wave components as 



Af = A'' 



Af 



(6) 



From Fig. 2(a), one can see that the d-wave order- 
parameter is suppressed within an about 5a radius 
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around the impurity site, with a being the lattice con- 
stant. Simultaneously, as shown in Fig. 2(b), an s-wave 
component is slightly induced with the magnitude of less 
than 1% relative to the bulk d-wave component. We note 
that the calculation with fixed ^ = ^q ( value for the uni- 
form case ) overestimates the induction of the s-wave 
compcpent up to 8%._pifferent from the cases with a 
vortexES' and surfaces |13^ the d- and s-wave components 
obtained here are real. The induced s-wave component 
has the following iJiaracteristics similar to the results of 
the previous workllj' ; (i) it has the same magnitude and 
opposite sign with tt/2 rotation centered on the impu- 
rity site, (ii) it is not induced on the sites in the diagonal 
direction from the impurity. These characteristics origi- 
nate from the sign-changing property of the (i2.2„j^2-wave 
order parameter, and the symmetry around the impurity. 



(a) 



rf-wave component 




(b) 



5-wave component 




Fig. 2. Spatial dependence of the relative amplitudes of (a) d- 
wave and (b) s-wave components of the order parameters in the 
unit cell illustrated in Fig. 1. Note the different scales for (a) 
and (b) . The order parameters are normalized by the bulk value 
of the d-wave component. The doping rate is 5 = 0.05 and 
J/t = 0.2. 



Next, we calculate the local density of states (LDOS) 
defined as 
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where i represents a site, k is the Bloch wave number 
and a is the eigenstate number. Though eq. (7) is calcu- 



lated within the Gutzwilkr approximation, we can ob- 
tain qualitative features.llS' Figure 3 (a) shows the LDOS 
obtained on the site- A, -B and -C illustrated in Fig 1. 
Here, we take Nc — 240 (15 x 16). On the sites located on 
the edge of the unit cell (site-C), we actually find the typ- 
ical DOS for d-wave superconductivity. Therefore, the 
impurity is considered to be well-isolated in the present 
situation; N^ = 27, S = 0.05. On the nearest-neighbor 
site of the impurity (site-^), the LDOS detects the zero- 
energy states (ZES) as the zero-energy peak (ZEP). This 
is due to the fact that, in the unitary scattering limit, 
there always exist some impurity-scattering processes in 
which the quasiparticles feel the sign-change of the pair 
potential. The asymmetric shape of the ZEE is due to 
the breaking of the electron-hole symmetry .c3 

Different from the case with vortex or surface, the 
splitting of the ZEP does not appear within the present 
calculation since the induced s-wave component is real 
and rather small. On the 3rd-neighbor site of the impu- 
rity (site-i3), the ZEP cannot be seen and the weight in 
the low energy region is reduced compared with the bulk 
d-wave case (site-C). This is due to the spatial variation 
of the ZES. In order to see it, the DOS with a fixed en- 
ergy _E = is plotted as a function of sites around the 
impurity in Fig. 3(b). This figure represents the spatial 
oscillating behavior of ZES, which can be regarded as the 
Friedel oscillation. This is the reason for the absence of 
the ZEP on site-B. Furthermore, Fig. 3(b) shows that 
there are no i-tails along the diagonal direction, and that 
the ZES is approximately localized around the impurity. 
This behavior is consistent with the results in refs. 6 and 
13 but is contrary to the results in refs. 9 and 11. Our 
preliminary results show that whether the ZES localizes 
or not depends on thainteraction strength and the shape 
of the Fermi surfaceO 

An interesting question is how the ZES behaves in the 
overdoped region. To examine the doping dependence 
of the ZES, attention should be given to the size of the 
unit cell since the coherence length becomes longer with 
increasing S. Even in the case oi 6 — 0.05, we find that 
small-size calculations {Nl < 23) give wrong results, i.e., 
they sometimes give a splitting of the ZEP. If we in- 
crease the size of the unit cell, the width of the splitting 
shrinks roughly in proportion to l/iVi, and the split- 
ting vanishes when N^ > 25. This size-dependence is 
due to the intercell overlapping effect of ZES and is not 
intrinsic to the single impurity case. Thus the appear- 
ance of the splitting of ZEP found by Tanaka et alr3 is 
due to the small cluster size they used {Nl = 16 ^ 18). 
Here, let us estimate the size of the unit cell required in 
the case oi 5 = 0.1. The obtained energy gap without 
any impurity is A = 0.27t at S = 0.05 and A = 0.22t 
at S — 0.1. The Fermi energy ep is estimated from 
ep ^ 4:{gtt+lgsJ\^Q\). We obtain ep ^ 0.79t at S = 0.05 
and ep ~ 1.1 It at 5 ~ 0.1. Then the ratio of the co- 
herence length at S — 0.05 to (5 = 0.1 is found to be 
{^^\s=o.i) / {^^\s=o.a5) ^ 1-5. Thus, we can estimate 
the size of the unit cell required for 5 = 0.1 as Np ^ 37. 
Indeed, our preliminary calculation with Np — 37 shows 
an almost single peak near the zero-energy in the LDOS. 
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Fig. 3. (a) The local density of states at the three sites A, B 
and C (see Fig.l), for & = 0.05 and J/t = 0.2. The size of the 
unit cell is 27 X 27 sites. Ni{E) has a zero-energy peak on the 
nearest-neighbor sites of the impurity site (site- A). We can see 
that Ni(E) reproduces the shape of the uniform d^2_„2-wave 
state on the sites located on the edge of the unit cell (site-C). 
(b) Spatial variations of the local density of states at i? = 
plotted over 17 X 17 sites around the impurity. 



Although this suggests the similar results in (5 = 0.1, the 
doping dependence of the ZES fromjjjider- to overdoped 
region is to be published elsewhere.E^ 

In this letter, we have shown for the first time that the 
ZES gives ZEP in the LDOS around a nonmagnetic im- 
purity in the (ia;2_j,2-wave superconducting state of the 
underdoped t-J model. In the scanning tunneling spec- 
troscopy (STS) experiment where the direct observation 
of the LDOS is possible with atomic-scale spatial resolu- 
tion, the ZBCP, which corresponds to the ZEP obtained 
here, is expected to be ocobed just as in the case of a well- 
oriented (110) surface.llZP Furthermore we found that 
the ZES is approximately localized (radius 5a) around 
the impurity. The induced s-wave component is real but 
small and thus there is no splitting of the ZEP. 
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